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We obtain strangeness production as function of centrality in a statis- 
tical hadronization model analysis of all experimental hadron production 
data in Au-Au collisions at ^/SNN = 62.4 GeV. Our analysis describes suc- 
cessfully the yield of strange and multi-strange hadrons recently published. 
We explore condition of hadronization as a function of centrality and find 
universality for the case of chemical non-equilibrium in the hadron phase 
space corresponding to quark-gluon plasma (QGP) in chemical equilibrium. 

PACS numbers: 24.10.Pa, 12.38.Mh, 25.75.-q, 13.60.Rj 

1. Experimental Data 

The statistical hadronization model (SHM) is known to successfully de- 
scribe hadron yields in a variety of experiments at different energies [I1I2]- 
Recently, new results on the yields of strange and multi-strange hadrons 
at ^/snn = 62.4 GeV have been published [,3j. This data together with 
previously published non-strange particle yields [1] and (f) meson measure- 
ment [5] provides a new opportunity to reevaluate hadronization conditions 
in ^/sNN = 62.4 GeV collisions. 

Heavier multistrange particles are more rarely produced and therefore 
to assure that statistical errors are small, they are often presented in wider 
centrality bins than non-strange particles. In order to study all particle 
yields as a function of centrality, we have to account for the different cen- 
trality binning. We associate each bin with average number of participants 
-^part- Then, we fit hadron yield as a function of A^'part using a power law, 

/(iVpart) = a-iVp^art + c, (1) 



(1) 
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where a, b and c are free parameters fitted and shown in table [TJ We see 
as expected that particles are always more numerous than their respective 
anti-particles. However for the fi, we had only three data points and three 
free parameters of the interpolating function from Eq. ([T|). Therefore, we 
fixed c = and fitted the other two. This fit, with on degree of freedom, 
gave us a qualitatively similar function as for the other particles. 

Table 1. Particle yields power fit parameters (as defined by Eq. ((TJ) used to describe 
particle yields as function of centrality. For ft, c=0 is assumed 
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2. SHM fit including multistrange hadrons 

We performed a cross check of SHARE [Ij, the thermal model imple- 
mentation we used, with three other codes currently in use by other groups 
using the compiled data in [2j. We tested our model on the data set from 
Au-Au collision at y/sNN = 200 GeV from STAR and obtained both yields 
and thermal parameters within 5% from the values in p]. 

Rather than hadron ratios studied in [3], we decided to fit hadron yields, 
which require a common volume dV/dy of the source. We complemented the 
data set with identified hadrons (vr^ , , ) from [4] and (f) yield from [5] . 
Furthermore, we require the conservation of charge per baryon {Q — Q)/{B — 
B) = 0.39ib0.01 and strangeness (s — s)/(s + s) = 0±0.05 with a few percent 
error to account for the detector acceptance and efficiency. 
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Fig. 1. Left panel: Particle yields as a function of number of participants with re- 
spective interpolating functions. Right panel: Top panel depicts chemical freeze- 
out temperature T as a function of the number of participants Npart ■ Red triangles 
are a copy of the results in I3_ , blue circles describe our fit of the same data set in 
semi-equilibrium model with 7g = 1, 7s 7^ 1 and purple squares are non-equilibrium 
model "fi l,i — q, s. Bottom panel shows the strangeness phase space occupancy 
compared to the light quark phase space occupancy ratio ^s/lq for different models. 
Open purple squares show the values of 7^ for the non-equilibrium model. 



As was also pointed out in [3j, the experimental particle yields are sub- 
ject to an important inconsistency: for all centralities, the yield of is 
smaller than both and (lowest of three corresponding lines in the 
left panel of the figured]). In general, at finite baryon density, we expect 
K+ <K° <K~. Interestingly the missing are associated with yields of 
A, A in that these three particle types had the largest contributions to the 
total of the fit. The choice to include or exclude them from the fit did 
not change the values of the SHM parameters we obtained. 

It is interesting to note that the three neutral particle types have the 
same V-type charged particle decay and that there is an overlap in kinematic 
ID of the decay channels. For this reason we fitted instead the yield of 
+ which does not have this ambiguity. Considering the smallness 
of the required shift from more numerous to A, A we also fitted the ratio 
A/A. These two data points replace the three raw yields and proved to be 
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consistent with all other particle yields. Assuming that other particles fix 
the source volume, this procedure incurs minimal loss of data information. 

We expected that our SHM fit with = 1 should reproduce results 
of [3]. And indeed, we obtain compatible chemical freeze-out temperature 
T (top right panel of figure [T]). However, we find that the strangeness phase 
space is strongly overpopulated, which is reflected by 7^ > 1 (bottom right 
panel of figure [1]) contrary to results published by STAR. After we found 
this discrepancy, we learned from the authors that in 7^ < 1 has been 
forced. This of course completely invalidates the conclusions of this work. 
This also resolves the myth that the yields of </) oc 7^ are inconsistent with 
SHM. 

As a next step, we considered the chemical non-equilibrium with ji 7^ 
1,1 = q,s. This approach allows that aside of the strangeness phase space 
also the light q = u, d quark phase space when observed on hadron side to be 
out of equilibrium. To preserve the number of degrees of freedom, we fitted 
the pressure P ~ 82 MeV /fm^ , which is the value of critical hadronization 
pressure found in SHM fits to other experiments [7j. 

The chemical non-equilibrium SHM fit results in an even more pro- 
nounced strangeness yield enhancement, 7^ — )• 1.5 at central rapidity, an 
over-population of the light quark phase space jg ~ 1.6, and a lower freeze- 
out temperature T ~ 142 MeV. Another important feature of this approach 
is that the parameter errors in the fit decreased dramatically compared to 
the semi-equilibrium model. In the right bottom panel of figure [H the cen- 
trality dependence of these parameters is shown and the errors on 7^ and 
7s for chemical non-equilibrium are within the size of the symbols. We note 
that the SHM fits at all centralities have confidence level above 60%. 

3. Hadronization conditions 

In the non-equilibrium SHM approach there appears to be universal 
hadronization condition as function of centrality, which can be seen in 
figure [2] top two panels on left. Apart from the constant pressure P ~ 
82MeV/fm^, obtained at other RHIC energies [7], which value was used in 
non-equilibrium fits and found consistent with the new data, the entropy 
density at hadronization stays at constant universal value a ~ 3.3 fm~^. 
Similarly, the energy density is constant at e ~ 0.48 GeV/fm'^. A universal 
value of statistical properties at hadronization as function of centrality [8] 
supports the notion of a new phase of matter, the quark-gluon plasma un- 
dergoing a sudden break up transition to hadrons which free stream, as the 
chemical non-equilibrium model assumes. On the other hand, in the semi- 
equilibrium hadronization model, both entropy and energy density grow 
nearly proportional to the number of participants. There is no physical 
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Fig. 2. Left panel: Going from top to bottom, centrality dependence of entropy 
density energy density (cr and e respectively) and strangeness per entropy ratio s/ S. 
Right panel: Going from top to bottom, centrality dependence of the strangeness 
rapidity density and energy per baryon. 



understanding of this behavior. 

We see, in the bottom panel on left in figure [21 that in SHM non- 
equilibrium fit, the ratio of strangeness s to entropy S approaches for most 
central collisions the QGP expectation s/S = 0.031. The relative yield of 
strangeness is systematically smaller in the SHM semi-equilibrium fit. On 
right in figure O we see at the top that strangeness density increases with 
centrality, in agreement with the expectation that the largest system can ap- 
proach the QGP yield. Moreover, we see that in the SHM non-equilibrium 
model, we always find a much greater strangeness yield compared to the 
SHM semi-equilibrium. The effect is quite large at small centrality, where 
the entropy contents differ most, see the left top panel. 

Another interesting result in non-equilibrium SHM is the constancy for 
-^part > 130 of hadronization energy per baryon E/B seen in bottom right 
of figure [21 Moreover, we find that the fraction of the available energy in 
the collision is universal among other systems and collision energies. We 
find, E/Braost central = 16.5 GcV ~ 0.25 X 62.4 GeV = jy/s^N. This result 
is the same as at the lower SPS energies. It says that the stopping of baryon 
number is 4 times more effective than the stopping of energy. Only at small 
centralities the stopping of energy and baryon number seems to converge to 
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a common value resulting in E/B^ost peripheral 2\/^^- 

4. Conclusions 

We revisited the statistical hadronization fit of particle yields from Au- 
Au collisions at ^Jsnn = 62.4 GeV recently updated with the yields of multi- 
strange hadrons. We obtained qualitatively different fit from the one pub- 
lished by the STAR collaboration [3J which was done with the constraints 
7g = 1 and 7^ < 1 and which fails to resolve yields of all multistrange 
particles. Our results show a very good confidence level at all centralities. 
Our results show that the hadron strangeness phase space is overpopulated 
7s > 1 from semi-central iVpart = 150 to most central collisions A'^part = 350. 
We allowed the light quark phase space, 7^ / 1, to be out of equilibrium 
as well and the new fit to the experimental data works for multistrange 
hadrons and has much smaller errors. However, given that were incon- 
sistent with the yields of charged hadrons, we did fit A -|- A -|- 2K^. We 
found that the non-equilibrium model provides physical consistency in that 
as a function of centrality, i.e., system size we find appropriate and universal 
behavior of bulk physical properties: entropy, energy, strangeness and the 
stopped energy per baryon [7]. 
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